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Abstract 

We consider a picture in which the transition from a big crunch to a big bang corresponds to the 



, collision of two empty orbifold planes approaching each other at a constant non-relativistic speed 



in a locally flat background space-time, a situation relevant to recently proposed cosmological 
■ models. We show that p-brane states which wind around the extra dimension propagate smoothly 

m . 

and unambiguously across the orbifold plane collision. In particular we calculate the quantum 



mechanical production of winding M2-branes extending from one orbifold to the other. We find 
that the resulting density is finite and that the resulting gravitational back-reaction is small. These 
winding states, which include the string theory graviton, can be propagated smoothly across the 



Qh! transition using a perturbative expansion in the membrane tension, an expansion which from 



the point of view of string theory is an expansion in inverse powers of a' . The conventional 
description of a crunch based on Einstein general relativity, involving Kasner or mixmaster behavior 
is misleading, we argue, because general relativity is only the leading order approximation to string 
theory in an expansion in positive powers of a'. In contrast, in the M theory setup we argue that 
interactions should be well-behaved because of the smooth evolution of the fields combined with the 
fact that the string coupling tends to zero at the crunch. The production of massive Kaluza-Klein 
states should also be exponentially suppressed for small collision speeds. We contrast this good 
behavior with that found in previous studies of strings in Lorentzian orbifolds. 
PACS number(s): 11.25.-w,04.50.+h, 98.80.Cq,98.80.-k 
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I. INTRODUCTION 

One of the greatest challenges faced by string and M theory is that of describing time- 
dependent singularities, such as occur in cosmology and in black holes. These singularities 
signal the catastrophic failure of general relativity at short distances, precisely the pathology 
that string theory is supposed to cure. Indeed string theory does succeed in removing the 
divergences present in perturbative quantum gravity about fiat spacetime. String theory 
is also known to tolerate singularities in certain static backgrounds such as orbifolds and 



2 



conifolds. However, studies within string theory thus far have been unable to shed much light 
on the far more interesting question of the physical resolution of time-dependent singularities. 

In thispaper we discuss M theory in one of the simplest possible time-dependent back- 
grounds (jj, |2|, a direct product of d — 1-dimensional flat Euclidean space R d ~ l with two 
dimensional compactified Milne space-time, hAc-, with line element 

-dt 2 + t 2 d6 2 . (1) 

The compactified coordinate 9 runs from to 9q. As t runs from — oo to +oo, the compact 
dimension shrinks away and reappears once more, with rapidity 9q. Analyticity in t suggests 
that this continuation is unique [3J. 

Away from t = 0, M.c is locally flat, as can be seen by changing to coordinates 
T = tcoshfl, Y = tsmh.9 in which JTJ) is just -dT 2 + dY 2 . Hence M c x is natu- 

rally a solution of any geometrical theory whose field equations are built from the curvature 
tensor. However, M. c x R d ~ l is nonetheless mathematically singular at t = because the 
metric degenerates when the compact dimension disappears. General relativity cannot make 
sense of this situation since there ceases to be enough Cauchy data to determine the future 
evolution of fields. In fact, the situation is worse than this: within general relativity, generic 
perturbations diverge as log|i| as one approaches the singularity |4j, signaling the break- 
down of perturbation theory and the approach to Kasner or mixmaster behavior, according 
to which the space-time curvature diverges as t~ 2 . Of course, this breakdown of general 
relativity presents a challenge: can M theory make sense of the singularity at t = 0? 

We are interested in what happens in the immediate vicinity of t — 0, when the compact 
dimension approaches, and becomes smaller than, the fundamental membrane tension scale. 
The key difference between M theory (or string theory) and local field theories such as general 
relativity is the existence of extended objects including those stretching across compactified 
dimensions. Such states become very light as the compact dimensions shrink below the 
fundamental scale. They are known to play a central role in resolving singularities for 
example in orbifolds and in topology-changing transitions |5|]. Therefore, it is very natural 
to ask what role such states play in big crunch/big bang space-times. 

In this paper we shall show that p-branes winding uniformly around the compact dimen- 
sion obey equations, obtained by canonical methods, which are completely regular at t = 0. 
These methods are naturally invariant under choices of worldvolume coordinates. Therefore 
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we claim that it is possible to unambiguously describe evolution of such states from t < to 
t > 0, through a cosmo logical singularity from the point of view of the low energy effective 
theory. Indeed the space-time we consider corresponds locally to one where two empty, 
flat, parallel orbifold planes collide, precisely the situation envisaged in recently proposed 
cosmological models. 

Hen. the captions we r epo rt are directly relevant to the erotic Q and cyclic 
universe |7( scenarios, in which passage through a singularity of this type is taken to represent 
the standard hot big bang. In particular, the equation of state during the dark energy and 
contracting phases causes the orbifold planes to be empty, flat, and parallel as they approach 
to within a string length 0,0]. This setup makes it natural to split the study of the collision 
of orbifold planes into a separate analysis of the winding modes, which become light near 
t = 0, and other modes that become heavy there. This strategy feeds directly into the 
considerations in this paper. 

The physical reason why winding states are well-behaved is easy to understand. The 
obvious problem with a space-time such as compactified Milne is the blue shifting effect felt 
by particles which can run around the compact dimension as it shrinks away. As we shall 
discuss in detail, winding states wrapping around the compact dimension do not feel any 
blue shifting effect because there is no physical motion along their length. Instead, as their 
length disappears, from the point of view of the noncompact dimensions, their effective mass 
or tension tends to zero but their energy and momentum remain finite. When such states 
are quantized the corresponding fields are well-behaved and the field equations are analytic 
at t — 0. In contrast, for bulk, non- winding states, the motion in the 9 direction is physical 
and it becomes singular as t tends to zero. In the quantum field theory of such states, this 
behavior results in logarithmic divergences of the fields near t = 0, even for the lowest modes 
of the field which are uniform in 9 i.e., the lowest Kaluza-Klein modes (see Section X and 
Appendix 4). 

We are specially interested in the case of M theory, considered as the theory of branes. 
As the compact dimension becomes small, the winding M2-branes we focus on are the 
lowest energy states of the theory, and describe a string theory in a certain time-dependent 
background. The most remarkable feature of this setup is that the string theory includes a 
graviton and, hence, describes perturbative gravity near t = 0. In this paper, we show these 
strings, when considered as winding M2-branes, follow smooth evolution (see Section VII) 
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across the singularity, even though the string frame metric degenerates there. Furthermore 
we show that this good behavior is only seen in a perturbation expansion in the membrane 
tension, corresponding from the string theory point of view to an expansion in inverse powers 
of a'. We argue that the two-dimensional nonlinear sigma model describing this situation 
is renormalizable in such an expansion. The good behavior of the relevant string theory 
contrasts sharply with the bad behavior of general relativity. There is no contradiction, 
however, because general relativity is only the first approximation to string theory in an 
expansion in positive powers of a'. Such an expansion is valid when t is much larger than 
the fundamental membrane scale, but it fails near the singularity where, as mentioned, 
the theory is regular in the opposite expansion. The logarithmic divergences of 

perturbations found using the Einstein equations are, thereby, seen to be due to the failure 
of the a' expansion, and not of M or string theory per se. 

When the M theory dimension is small, the modes of the theory are neatly partitioned 
into light ^-independent modes and heavy ^-dependent modes. The former set consists of 
winding membranes, which describe a string theory including perturbative gravity. This is 
the sector within which cosmological perturbations lie, and which will be our prime focus in 
this paper. 

The ^-dependent modes are likely to be harder to describe. The naive argument that 
these modes are problematic because they are blue shifted and, hence, infinitely amplified as 
t — > is suspect because it relies on conventional Einstein gravity, Here we argue that, close 
to the brane collision, Einstein gravity is a poor approximation and, instead, perturbative 
gravity is described by the non- singular winding sector. The latter does not exhibit blue 
shifting behavior near t = 0, so the naive argument does not apply. 

Witten has argued that the massive Kaluza-Klein modes of the eleven-dimensional 
theory map onto non-perturbative black hole states in the effective string theory. Even 
though these black hole states are likely to be hard to describe in detail, we will explain in 
Section II why their overall effect is likely to be small. First, in the cosmological scenarios 
of interest, the universe enters the regime where perturbative gravity is described by the 
winding modes (i.e., the branes are close) with a negligible density of Kaluza-Klein massive 
modes. This suppression is a result of the special equation of state in the contracting phase 



that precedes this regime 



Second, the density of black holes quantum produced due 



to the time-dependent background in the vicinity of the collision is likely to be negligible 
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because they are so massive and so large. 

For these reasons, we focus at present on the propagation of the perturbative gravity sec- 
tor near t = corresponding to the winding M2-brane states. In Section II, we introduce the 
compactified Milne background metric that describes the collision between orbifold planes 
in the big crunch/big bang transition. We also discuss the motivation for the initial condi- 
tions that will be assumed in this paper. The canonical Hamiltonian description of p-branes 
in curved space is given in Section III and applied to winding modes in the compactified 
Milne background in Section IV. Section V discusses the key difference in the Hamiltonian 
description between winding and bulk states that accounts for their different behavior near 
the big crunch/big bang transition. 

Section VI is the consideration of a toy model in which winding strings are produced as 
the branes collide. The winding modes are described semi-classically, and their quantum 
production at the bounce is computed. Section VII presents the analogous semi-classical 
description of winding M2-branes. Although we cannot solve the theory exactly, we show 
the eleven-dimensional theory is well behaved near t = and explain how the apparent 
singularity in the dimensionally-reduced string theory is resolved in the membrane picture. 
Then, Section VIII makes clear the difference between our calculation, an expansion in 
inverse powers of a', versus Einstein gravity, the leading term in an expansion in positive 
powers of a'. This argument is key to explaining why we think the transition is calculable 
even though it appears to be poorly behaved when described by Einstein gravity. Section 
IX, then, uses Euclidean instanton methods to study the quantum production of winding 
M2-branes (in analogy to the case of winding strings in Section VI) induced by passage 
through the singularity, obtaining finite and physically sensible results. In particular, the 
resultant density tends to zero as the speed of contraction of the compact dimension is 
reduced. We estimate the gravitational back- reaction and show it is small provided 80, the 
rapidity of contraction of the compact dimension, is small. 

In Section X, we comment on why our M theory setup is better behaved than the 
Lorentzian orbifold case Q| considered in some previous investigations of the big crunch/big 
bang transition. The fundamental problem with the latter case, we argue, is that perturba- 
tive gravity lies within the bulk sector and not the winding sector as far as the compactified 
Milne singularity is concerned. Therefore, it is susceptible to the blueshifting problem men- 
tioned above, rendering the string equations singular. 
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II. THE BACKGROUND BIG CRUNCH/BIG BANG SPACE-TIME 



The d + 1- dimensional space-time we consider is a direct product of d — 1-dimensional 
Euclidean space, R d ~ l , and a two-dimensional time-dependent space-time known as com- 
pactified Milne space-time, or M.q- The line element for Aic x R d ~ x is thus 

ds 2 = -dt 2 + t 2 d6 2 + dx 2 , 0<9<9 , -oo<t<oo, (2) 

where x are Euclidean coordinates on 9 parameterizes the compact dimension and t 

is the time. The compact dimension may either be a circle, in which case we identify 9 with 
9 -\- 9q, or a. Z2 orbifold in which case we identify 9 with 9 + 29q and further identify 9 with 
29q — 9. The fixed points 9 = and 9 = 9q are then interpreted as tensionless Z2-branes 
approaching at rapidity 9o, colliding at t = to re-emerge with the same relative rapidity. 

The orbifold reduction is the case of prime interest in the ekpyrotic/cyclic models, orig- 
inally motivated by the construction of heterotic M theory from eleven dimensional super- 
gravity 10, In these models, the boundary branes possess nonzero tension. However, 
the tension is a subdominant effect near t = and the brane collision is locally well-modeled 
by M c x R d ~ l (See Ref. Q) 

The line element (0) is of particular interest because it is locally flat and, hence, an exact 
solution not only of d+ 1-dimensional Einstein gravity but of any higher dimensional gravity 
theory whose field equations are constructed from curvature invariants with no cosmological 
constant. And even if a small cosmological constant were present, it would not have a large 
effect locally so that solutions with a similar local structure in the vicinity of the singularity 
would be expected to exist. 

Consider the description of (j2J) within d + 1-dimensional general relativity. When the 
compact dimension is small, ^-dependent states become massive and it makes sense to 
describe the system using a low energy effective field theory. This may be obtained by the 
well known procedure of dimensional reduction. The d + 1-dimensional line element (j2J) may 
be rewritten in terms of a d- dimensional Einstein frame metric, gffi, and a scalar field 0: 



ds 2 = e 2*V(<i-2)/(d-l) d 02 + e -2H^(d-2) { d-l) g{ d) dx » dx v_ (3) 

The numerical coefficients are chosen so that if one substitutes this metric into the d + 1- 
dimensional Einstein action and assumes that cf> and gjfj are both ^-independent, one obtains 
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(i-dimensional Einstein gravity with a canonically normalized massless, minimally coupled 
scalar field <fi. (Here we choose units in which the coefficient of the Ricci scalar in the d- 
dimensional Einstein action is ^. We have also ignored Kaluza-Klein vectors, which play no 
role in this argument and are in any case projected out in the orbifold reduction.) 

From the viewpoint of the low energy effective theory, the d + 1-dimensional space-time 
M.c x R d ~ l is reinterpreted as a (i-dimensional cosmological solution where t plays the role 
of the conformal time. Comparing (J2J) and (jUJ), the d-dimensional Einstein- frame metric 



9$ = a2r lfj,is with a oc \t\ l ^ d ~ 2 ^ and the scalar field <fi = J{d— l)/(d — 2)ln|t|. From this 
point of view t = is a space-like curvature singularity of the standard big bang type where 
the scalar field diverges, and passing through t = would seem to be impossible. However, 
by lifting to the higher dimensional viewpoint one sees that the situation is not really so 
bad. The line element (J2J is in fact static at all times in the noncompact directions x. So for 
example, matter localized on the branes would see no blue shifting effect as the singularity 
approach^ . As we discuss in detail in Section IV, winding states do not see a blue shifting 
effect either. 

In this paper, we consider an M theory picture with two empty, flat, parallel colliding 
orbifold planes and we are interested in the dynamics of the collision region from the point 
where the planes are roughly a string length apart. The assumed initial conditions are im- 
portant for two reasons. First, they correspond to the simple compactified Milne background 
discussed above. Second, as mentioned in the introduction, this initial condition means that 
the excitations neatly divide into light winding modes that are becoming massless and heavy 
Kaluza-Klein modes that are becoming massive and decoupling from the low-energy effective 
theory. 

What we want to show now is that initial conditions with negligible heavy Kaluza-Klein 
modes present are naturally produced in cosmological scenarios such a the cyclic model [?], 
This justifies our focus on the winding modes throughout the remainder of the paper. 
However, the argument is inessential to the rest of the paper and readers willing to accept 
the initial conditions without justification may wish to proceed straight away to the next 
Section. 

The cyclic model assumes a non-perturbative potential hat produces an attractive force 
between the orbifold planes. When the branes are far apart, perhaps 10 4 Planck lengths, the 
potential energy is positive and small, acting as the dark energy that causes the currently 
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observed accelerated expansion. In the dark energy dominated phase, the branes stretch by 
a factor of two in linear dimensions every 14 billion years or so, causing the branes to become 
flat, parallel and empty. In the low-energy effective theory, the total energy is dominated 
by the scalar field <fi whose value determines the distance between branes. As the planes 
draw together, the potential energy V(4>) of this field decreases and becomes increasingly 
negative until the expansion stops and a contracting phase begins. 

A key point is that this contracting phase is described by an attractor solution, which 
has an equation of state parameter w = P/p 3> 1. The energy density of the scalar field 
scales as 

P4> oc a -(^-i)(i+») (4) 

in this phase. This is a very rapid increase, causing the density in <fi to come to dominate 
over curvature, anisotropy, matter, or radiation 9]. We now show that d> comes to dominate 
over the massive Kaluza-Klein modes. The latter scale as 



Pkk oc a 



-(d-l) L -(d-l)/(d-2) ^ 



where L is the size of the extra dimension. The first factor is the familiar inverse volume 
scaling which all particles suffer. The second factor indicates the effective mass of the 
Kaluza-Klein modes. The d + 1-dimensional mass is L _1 , but this must be converted to 
a d— dimensional mass using the ratios of square roots of the 00 components of the d + 1- 
dimensional metric and the (i-dimensional metric. This correction produces the second factor 
in ©. 

From Q, we have L cx e 0\A d ~ 2 )/( d-1 ) . Now the key point is that this scales much more 
slowly with a than the potential V(4>) which scales as p^ in (J3J). Neglecting the scaling with 
L, the density of massive Kaluza-Klein modes scales as as pi . The final suppression 
of the density of massive modes relative to the density in <p is therefore ~ (yi/Vf) w ^ l+w " > 
where Vi and Vf are the magnitudes of the scalar potential when the w ^> 1 phase begins 
and ends. For large w, which we need in order to obtain scale-invariant perturbations, this 
is an exponentially large factor j?l 0| . 

The massive Kaluza-Klein modes are, hence, exponentially diluted when the w ^> 1 phase 
ends and the Milne phase begins. During the Milne phase, the scalar field is massless and 
has an equation of state w = 1, so p^ scales as a" 2 ^ 1 ^ as the distance between the orbifold 
planes shrinks to zero. In this regime, the Kaluza-Klein massive mode density scales in 



precisely the same way. Therefore, their density remains an exponentially small fraction of 
the total density right up to collision: meaning from the string theory point of view that 
the black hole states remain exponentially rare. 

So we need only worry about black holes produced in the vicinity of the brane collision 
itself. From the point of view of the higher dimensional theory, the oscillation frequency 
of the masssive Kaluza-Klein modes u ~ 1 6*0^ I 1 changes adiabatically, lo/uj 2 ~ 8® <C 1 for 
small #o, all the way to t — 0. Therefore, one expects little particle production before or after 
t = 0. From the dimensionally reduced point of view, the mass of the string theory black 
holes is larger than the Hubble constant ~ t -1 , by the same factor 9q . From either analysis, 
production of such states should be suppressed by a factor e^ 1 / 6 * , making it negligible for 
small 9 . 

In sum, for the cosmological models of interest, the Kaluza-Klein modes are exponentially 
rare when the Milne phase begins, and, since their mass increases as the collision approaches, 
they should not be generated by the orbifold plane motion. (They effectively decouple from 
the low energy effective theory.) Hence, all the properties we want at the outset of our 
calculation here are naturally achieved by the contracting phase with w > 1, as occurs in 
some current cosmological models ^] . 

III. GENERAL HAMILTONIAN FOR p-BRANES IN CURVED SPACE 

The classical and quantum dynamics of p-branes may be treated using canonical methods, 
indeed p-branes provide an application par excellence of Dirac's general method. As Dirac 
himself emphasized |l9j, one of the advantages of the canonical approach is that it allows a 
completely general choice of gauge. In contrast, gauge fixed methods tie one to a choice of 
gauge before it is apparent whether that gauge is or isn't a good choice. In the situation of 
interest here, where the background space-time is singular, the question of gauge choice is 
especially delicate. Hence, the canonical approach is preferable. 

In this Section we provide an overview of the main results. The technical details are 
relegated to Appendix 1. Our starting point is the Polyakov action for a p-brane described 
by embedding coordinates in a a background space-time with metric g^: 



where n P is a mass per unit p-volume. The p-brane worldvolume has coordinates a a , where 
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cr° = r is the time and a 1 , i = l..p are the spatial coordinates. 

Variation of the action with respect to 7 a/ 3 yields the constraint that for p ^ 1, ^y a p 
equals the induced metric d a x ,1 dpx l 'g liv whereas for p = 1 j a p is conformal to the induced 
metric. Substituting these results back into the action one obtains the Nambu action for 
the embedding coordinates x M (cr Q ) i.e., —\i p times the induced p-brane world volume. We 
shall go back and forth between the Polyakov and Nambu forms in this paper. The former 
is preferable for quantization but the latter is still useful for discussing classical solutions. 

The simplest case of (0) is p = 0, a 0-brane or massive particle. Writing 7 o = — e 2 with 
e the 'einbein', one obtains 

1 r 

where we have set /Iq = m and the dot above a variable indicates a derivative with respect 
to r. Variation with respect to e yields the constraint e 2 = —x^x u g^ v . The canonical 
momentum is p M = rag n U x p 'e~ l and the constraint implies the familiar mass shell condition 

g^PfiPu = tv? . 

The canonical treatment for general p is explained in Appendix 1. The main result is 
that a p-brane obeys p + 1 constraints, reading 

C = K^gT + //jDet^x^) w 0, Q = w 0, (8) 

where '~ 0' means 'weakly zero' in sense of the Dirac canonical procedure (see Appendix 
I). Here the brane embedding coordinates are x^ and their conjugate momentum densities 
are tt^. The spatial worldvolume coordinates are a 1 , i = l,...,p, and the corresponding 
partial derivatives are denoted x^. The quantity x^x^g^ is the induced spatial metric on 
the p-brane. In Appendix 2 we calculate the Poisson bracket algebra of the constraints (JHJ), 
showing that the algebra closes and hence the constraints are all first class. The constraints 
(JHJ) are invariant under worldvolume coordinate transformations. 

The Hamiltonian giving the most general evolution in worldvolume time r is then given 

by 

H = J <Fa (^AC + A*C^ , (9) 

with C and C, given in (JHJ). The functions A and A 1 are completely arbitrary, reflecting 
the arbitrariness in the choice of worldvolume time and space coordinates. All coordinate 
choices related by nonsingular coordinate transformations give equivalent physical results. 
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For p = 0, anything with a spatial index i can be ignored, except the determinant in 
(JHJ) which is replaced by unity. The first constraint is then the usual mass shell condition, 
and the Hamiltonian is an arbitrary function of r times the constraint. The case of p = 1, 
i.e., a string, in Minkowski space-time, = rj^ is also simple and familiar. In this 
case, the constraints and the Hamiltonian JS} are quadratic. The resulting equations of 
motion are linear and hence exactly solvable. The constraints (jHJ) amount to the usual 
Virasoro conditions. In general, the p + 1 constraints (jHJ) together with the p + 1 free choices 
of gauge functions A and reduce the number of physical coordinates and momenta to 
2(d + 1) — 2(p + 1) = 2(d — p), the correct number of transverse degrees of freedom for a 
p-brane in d spatial dimensions. 

IV. WINDING p-BRANES IN Mc x R^ 1 

In this paper, we shall study the dynamics of branes which wind around the compact 
dimension in M.q x R , the line element for which is given in (J2J). This space-time possesses 
an isometry 9 — > #+constant, so one can consistently truncate the theory to consider p- 
branes which wind uniformly around the 9 direction. Such configurations may be described 
by identifying one of the p-brane spatial coordinates (the p'th spatial coordinate, a p say) 
with 9 and to simultaneously insist that that dpX^ 1 = d p ir^ = 0. 

Through Hamilton's equations, the constraint 9 = a p implies that ng = 0. This suggests 
that we can set 9 = a p and irg = and, hence, dimensionally reduce the p-brane to a (p — 1)- 
brane. Detailed confirmation that this is indeed consistent proceeds as follows. We compute 
the Poisson brackets between all the constraints C, Ci, 9 — a v and rcg. Following the Dirac 
procedure, we then attempt to build a maximal set of first class constraints. The constraints 
C and Ci commute with each other, for all a 1 , but not with 9 — a p and ttq. The solution 
is to remove all the irg and # p terms from C and Cg by adding terms involving ng and 
9 )P — 1 = (9 — cr p ), p . The new C and C, are now first class since they have weakly vanishing 
Poisson brackets with all the constraints, and the remaining second class constraints are 
9 — a p and TTg. For these, construction of the Dirac bracket is trivial and it amounts simply 
to canceling the 9 and Tig derivatives from the Poisson bracket. The conclusion is that we 
can indeed consistently set 9 = a p and ng = 0. We shall see in the following section that 
eliminating 9 and 7Tg in this way results directly in the good behavior of the winding modes 
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as t — > 0, in contrast with the bad behavior of bulk modes. 

The surviving first class constraints for winding p-branes are those obtained by substi- 
tuting 9 = a p and iig = into the p-brane constraints (jHI), namely 

C = tt^tT + /i^ 2 Det(z%.7^) ^ 0; C { = xfa « 0, (10) 

where i and j now run from 1 to p — 1 and fi and z/ from to d. The t 2 term comes from 
the 99 component of the M. c x R d ~ x background metric (0). We have also re-defined the 
momentum density ir^ for the p — 1-brane to be 9q times the momentum density for the p- 
brane so that the new Poisson brackets are correctly normalized to give a p — 1-dimensional 
delta function. The Hamiltonian is again given by the form © with the integral taken over 
the remaining p — 1 spatial coordinates. 

For p — 1, the reduced string is a <i-dimensional particle, ir^ is now the momentum p M 
and the determinant appearing in (|1U|) should be interpreted as unity. The second constraint 
is trivial since there are no remaining spatial directions. The general Hamiltonian reads: 

H = A(t) + »W) , (11) 

where fi, v run from to d — 1 and A{r) is an arbitrary function of r. We shall study the 
quantum field theory for this Hamiltonian in Section VI. 

Comparing (fTUj) with (JHJ) , we see that a p-brane which winds around the compact dimen- 
sion in Aic x R d ~ x behaves like a p — 1-brane in Minkowski spacetime with a time-dependent 
effective tension fi p 9o\t\, i.e., the p-brane tension times the size of the compact dimension, 
9 \t\. 

V. WINDING STATES VERSUS BULK STATES 

We have discussed in detail how in the canonical treatment the coordinate 9 and conjugate 
momentum density ng may be eliminated for p-branes winding uniformly around the compact 
dimension. This is physically reasonable, since motion of a winding p-brane along its own 
length (i.e. along 9) is meaningless. This is a crucial difference from bulk states. Whereas 
the metric on the space of coordinates for bulk states includes the t 2 d9 2 term, the metric on 
the space of coordinates for winding states does not. As we discuss in detail in Appendix 
4, when we quantize the system the square root of the determinant of the metric on the 
space of coordinates appears in the quantum field Hamiltonian. For bulk modes the metric 
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on the space of coordinates inherits the singular behavior of the background metric (J2J), 
degenerating at t = so that causing the field equations to become singular at t — even for 
^-independent field modes (see Section X). Conversely, for winding modes the Hamiltonian 
operator is regular at t — 0. 

The metric on the space of coordinates is defined by the kinetic energy term in the action: 
if the action reads S = ~ / drgjjx 1 ^ + . . ., where x 1 are the coordinates, then gu is the 
metric on the space of coordinates. The sum over / includes integration over a in our case. 
This superspace metric is needed for quantizing the theory, for example in the coordinate 
representation one needs an inner product on Hilbert space and this involves integration over 
coordinates. The determinant of gu is needed in order to define this integral (see Appendix 
4). 

The simplest way to identify the physical degrees of freedom is to choose a gauge, for 
example A=constant, = 0. For bulk particles, we can then read off the metric on 
coordinate space from the action (J2J) - in this case it is simply the background metric itself. 

We have already derived the Hamiltonian for winding states, and showed how through the 
use of Dirac brackets the 9 coordinate may be discarded. If we choose the gauge A = 1, Ai = 
in the Hamiltonian Q with constraints given in ([10)1 . we can construct the corresponding 
gauge-fixed action: 

S gf = J dr<F- l a l - - /iJ_^ 2 t 2 Det(x>^^)) (12) 

where fi, v run over to d — 1 and i,j run from 1 to p — 1. One may check that the 
classical equations following from the action ()12j) are the correct Lagrangian equations for 
the p-brane in a certain worldvolume coordinate system and that these equations preserve 
the constraints (jTUj) (see Appendix 3). 

The metric on the space of coordinates may be inferred from the kinetic term in (|12j). 
and it is just the Minkowski metric. In contrast, as discussed, the metric on the space of 
coordinates for bulk states involves the full background metric (J2J) which degenerates at 
t = 0. The difference means that whereas the quantum fields describing winding states 
are regular in the neighborhood of t = 0, those describing bulk states exhibit logarithmic 
divergences. In the penultimate section of this paper we argue that these divergences are 
plausibly the origin of the bad perturbative behavior displayed by strings and particles 
propagating on Lorentzian orbifolds, behavior we do not expect to be exhibited in M2-brane 
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winding states in M theory. 



VI. TOY MODEL: WINDING STRINGS IN M c x i? d 1 

Before approaching the problem of quantizing winding membranes, we start with a toy 
model consisting of winding string states propagating in M. c x R 4 ^ 1 ■ This problem has also 



been considered by others |l6l . 1 171 ] and in more detail than we shall do here. They point 
out and exploit interesting analogies with open strings in an electric field. Our focus will be 
somewhat different and will serve mainly as a warmup for case of winding M2-branes which 
we are more interested in. 

Strings winding uniformly around the compact 9 dimension in (J2J) appear as particles 
from the d— dimensional point of view. To study the classical behavior of these particles, it 
is convenient to start from the Nambu action for the string, 

S = -fiJ <P<T<J-T>et(d a xi*d f ,x''g lu ,) , (13) 

where \x is the string tension (to avoid clutter we set \i\ = \i for the remainder of this 
section). The string worldsheet coordinates are a a = (r, a). 

For the winding states we consider, we can set 9 = er, so < a < 6 . We insist that the 
other space-time coordinates of the string x M = (t, x) do not depend on a. It is convenient 
also to choose the gauge t = r, in which the action (JT3J) reduces to 



S = -1260 J dt\t\yi -x 2 , (14) 

in which t is now the time, not a coordinate. This is the usual square root action for a 
relativistic particle, but with a time-dependent mass fi6o\t\- The canonical momentum is 



p = fi9o\t\x/yl — x 2 and the classical Hamiltonian generating evolution in the time t is 



H = yp 2 + (n9ot) 2 . This is regular at t = 0, indicating that the classical equations should 
be regular there. 

Due to translation invariance, the canonical momentum p is a constant of the motion. 
Using this, one obtains the general solution 

P 

x = x -\ — — sinh -1 (fj,8 t/\p\) , — oo < t < oo, (15) 

according to which the particle moves smoothly through the singularity. At early and late 
times the large mass slows the motion to a crawl. However, at t — the particle's mass 
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disappears and it instantaneously reaches the speed of light. The key point for us is that 
these winding states have completely unambiguous evolution across t — 0, even though the 
background metric (|2J) is singular there. 

Now we turn to quantizing the theory, as a warmup for the membranes we shall consider 
in the next section. The relevant classical Hamiltonian was given in it describes a 

point particle with a mass ^9o\t\. In a general background space-time, ordering ambiguities 
appear, which are reviewed in Appendix 4. However, in the case at hand, there are no 
such ambiguities. The metric on the space of coordinates is the Minkowski metric r]^ u . The 
standard expression for the momentum operator p v = —ih(d/dx u ), and the Hamiltonian H 
given in are clearly hermitian under integration over coordinate space, / d d x. Finally, 
the background curvature R vanishes for our background so there is no curvature term 
ambiguity either. 

Quantization now proceeds by setting p^ = —id^ (we use units in which h is unity) in 
the Hamiltonian constraint fTTf which is now an operator acting on the quantum field 4>. 
Fourier transforming with respect to x, we obtain 

= - (p 2 + (fx9ot) 2 ) 0, (16) 

i.e., the Klein Gordon equation for a particle with a mass ^9 \t\. 

Equation (fTTfl) is the parabolic cylinder equation. Its detailed properties are dis- 
cussed in Ref. [20], whose notation we follow. We write the time-dependent frequency 
as uj = \jp 2 + (fi9 t) 2 . At large times fi9 \t\ ^> \p\, u is slowly varying: ujjuj 1 <C 1 so 
all modes follow WKB evolution. The general solution behaves as a linear combination of 
u~hxp(±i J udt) « t - 5exp(±i~ (^9 t 2 + (p 2 lnt)/(/i6» )))). 

For large momentum, p 2 ^> [i9o, the WKB approximation remains valid for all time 
since uj/uj 2 is never large. In the WKB approximation there is no mode mixing and no 
particle production. Therefore for large momentum one expects little particle production. 
Departures from WKB are nonperturbative in uj/u 2 , as explicit calculation verifies, the 
result scaling as ~ exp(— \u 2 /u\ max ) ~ exp(—p 2 /fi9 ), at large p 2 . 

The parabolic cylinder functions which behave as positive and negative frequency modes 
at large times are denoted E(a, x) and E*(a, x), where x = \/2fi9 t and a = —p 2 / (2/i#o). For 
positive x they behave respectively as x~5 exp(+z(lx 2 — alna;) and x~^ exp(— i{\x 2 — alnx). 
Both E(a,x) and E*(a,x) are analytic at x — 0. They are uniquely continued to negative 
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values through the relation 

E(a, -x) = -ie™E{a, x) + Wl + e 2 ™E*(a, x). (17) 

For t < 0, E(a, —y^JiUot) is the positive frequency incoming mode. As we extend t to 
positive values, (|17|) yields the outgoing solution consisting of a linear combination of the 
positive frequency solution E* and the negative frequency solution E. The Bogoliubov 
coefficient |2l| (3 for modes of momentum p is read off from ()17j) : 

p = -ie™ = -i e -hP 2 l^\ (18) 

The result is exponentially suppressed at large p, hence, the total number of particles per 
unit volume created by passage through the singularity is 

which is finite and tends to zero as the rapidity of the brane collision is diminished. 

It is interesting to ask what happens if we attempt to attach the t < half of M.c x -R , 
(J2J) with rapidity parameter 6™, to the upper half of A4c x R d_1 with a different rapidity 
parameter 0Q Ut . After all, the field equations for general relativity break down at t = and, 
hence, there is insufficient Cauchy data to uniquely determine the solution to the future. 
Hence, it might seem that we have the freedom to attach a future compactified Milne with 
any parameter #q u ', since this would still be locally flat away from t = and, hence, a 
legitimate string theory background. However it is quickly seen that this is not allowed. 
By matching the field <fi and its first time derivative d t <p across t = 0, we can determine 
the particle production in this case. We find that due to the jump in 9 , the Bogoliubov 
coefficient (3 behaves like ((0™)% — (0Q lt )^) / \J 0Q n #o"*, at large momentum, independent of 
p. This implies divergent particle production, and indicates that to lowest order one only 
obtains sensible results in the analytically- continued background, which has 6q u1 = 6™. We 
conclude that to retain a physically sensible theory we must have 9q u1: = 9 l n , at least at 
lowest order in the interactions. 

VII. DYNAMICS OF WINDING M2-BRANES 

We now turn to the more complicated but far more interesting case of winding membranes 
in M.q x R d ~ l - We have in mind eleven-dimensional M theory, where the eleventh, M theory 
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dimension shrinks away to a point. When this dimension is small but static, well known 



arguments |22j indicate that M theory should tend to a string theory: type IIA for circle 
compactification, heterotic string theory for orbifold compactification. It is precisely the 
winding membrane states we are considering which map onto the string theory states as 
the M theory dimension becomes small. What makes this case specially interesting is that 
the string theory states include the graviton and the dilaton. Hence, by describing string 
propagation across t = we are describing the propagation of perturbative gravity across a 
singularity, which as explained in Section II is an FRW cosmological singularity from the d 
dimensional point of view. 

A winding 2-brane is a string from the <i-dimensional point of view. As explained in 
section III, the Hamiltonian for such strings may be expressed in a general gauge as 

H = I da ( \ + V^'^') + ^ V'vv) , (20) 

where A and A 1 are arbitrary functions of the worldsheet coordinates a and r, and prime 
represents the derivative with respect to a. Here //, v run over 0, l,..,d— 1, and primes 
denote derivatives with respect to a. The Hamiltonian is supplemented by the following 
first class constraints 

WvtT + ^ 2 2 o t\ux^x v/ « 0; « 0, (21) 

which ensure the Hamiltonian is weakly zero. The latter constraint is the familiar require- 
ment that the momentum density is normal to the string. 

The key point for us is that the Hamiltonian and the constraints are regular in the 
neighborhood of t = 0, implying that for a generic class of worldsheet coordinates the 
solutions of the equations of motion are regular there. 

As we did with particles, it is instructive to examine the classical theory from the point of 
view of the Nambu action. We may directly infer the classical action for winding membranes 
on M.c x R d ~ x by setting 9 = a 2 and d 2 t = d 2 x = 0. The Nambu action for the 2-brane 
then becomes 

S = -yMo / dVltl^-Det^x^x"^), (22) 
where a a = (r, a) and /i runs over 0, ..d — 1. This is precisely the action for a string (|13[) 



in a time- dependent background g^ u = 9 
corresponding to M. c x R d ~ x in M theory 



r]^, the appropriate string- frame background 



!]■ 
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As a prelude to quantization, let us discuss the classical evolution of winding membranes 
across t = 0. We can pick worldsheet coordinates so that x° = t = t, and x-x' — 0. In this 
gauge the Nambu action is: 

S = -fx 2 J dtda\t\\x'\ \fl - £ 2 , (23) 
and the classical equations are 

d t (ef) = d a {jd^ , d t e = t^, (24) 

where 

e = \t\^Jx' 2 /(l-S 2 ) (25) 

the energy density along the string is fi2t- It may be checked that the solutions of these 
equations are regular at t = 0, with the energy and momentum density being finite there. 
The local speed of the string hits unity at t — 0, but there is no ambiguity in the resulting 
solutions. The data at t — consists of the string coordinates x(o) and the momentum 
density vr(cr) = fi2tx, which must be normal to x'(a) but is otherwise arbitrary. This is the 
same amount of initial data as that pertaining at any other time. 

Equations ()25|) describe strings in a d dimensional flat FRW cosmological background 
with ds 2 = aft) 2 !]^ and scale factor a{t) oc . The usual cosmological intuition is helpful in 
understanding the string evolution. The comparison one must make is between the curvature 
scale on the string and the comoving Hubble radius \t\. When \t\ is larger than the comoving 
curvature scale on the string, the string oscillates as in flat space-time, with fixed proper 
amplitude and frequency. However, when \t\ falls below the string curvature scale, the 
string is 'frozen' in comoving coordinates. This point of view is useful in understanding the 
qualitative behavior we shall discuss in the next section. 

There is one final point we wish to emphasize. In Section II we discussed the general 
Hamiltonian for a p-brane in curved space. As we have just seen, a winding M2-brane on 
M.c x R d ~ l has the same action as a string in the background g^ u = \t\rjfu,. We could have 
considered this case directly using the methods of Section II. The constraints (jEJ) for this 
case read: 

(#o|t|ry%^ + ^\Ht\i]^'x v ' w 0; x^ w 0, (26) 

and the Hamiltonian would involve an arbitrary linear combination of the two. These 
constraints may be compared with those coming directly from our analysis of winding 2- 
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branes, given by (|10j) with p = 2. Only the first constraint differs, and only by multiplication 
by 6o\t\. For all nonzero 9o\t\, the difference is insignificant: the constraints are equivalent. 
Multiplication of the Hamiltonian by any function of the canonical variables merely amounts 
to a re-definition of worldsheet time. This is just as it should be: dimensionally reduced 
membranes are strings. However, the membrane viewpoint is superior in one respect, namely 
that the background metric (J2J) is non-singular in the physical coordinates (which do not 
include 9). That is why the membrane Hamiltonian constraint is nonsingular for these states. 
The membrane viewpoint tells us we should multiply the string Hamiltonian by 9o\t\ in order 
to obtain a string theory which is regular at t = 0. Without knowing about membranes, 
the naive reaction might have been to discard the string theory on the basis that the string 
frame metric is singular there. 

VIII. EINSTEIN GRAVITY VERSUS AN EXPANSION IN I /a! 

We are interested in the behavior of M theory, considered as a theory of M2-branes, in 
the vicinity of t = 0. The first question is, for what range of \t\ is the string description 
valid? The effective string tension \i\ is given in terms of the M2-brane tension \i2 by 

An = yMoN- (27) 

i 

The mass scale of stringy excitations is the string scale fif . For the stringy description to be 
valid, this scale must be smaller than the mass of Kaluza-Klein excitations, (# |£|) _1 - This 
condition reads 

\t\ < /i 2 %\ (28) 

The second question is: for what range of \t\ may the string theory be approximated by d- 
dimensional Einstein-dilaton gravity? Recall that the string frame metric is g^ v = a 2 r]^ v = 

\t\V)Jivi so the curvature scale is a/a ~ The approximation holds when the curvature 

i 

scale is smaller than the string scale p,\ , which implies 

\t\>^k~K (29) 

We conclude that for small 9q there are three relevant regimes. In units where the membrane 
tension /i 2 is unity, they are as follows. For \t\ > 9Q 1 the description of M theory in terms of 
eleven dimensional Einstein gravity should hold. As \t\ falls below 9q 1 , the size of the extra 
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dimension falls below the membrane scale and we go over to the ten dimensional string 

theory description. In the cosmo logical scenarios of interest, the incoming state is very 

smooth ^,0], hence this state should be well-described by ten dimensional Einstein-dilaton 

gravity, which of course agrees with the eleven dimensional Einstein theory reduced in the 

_ i 

Kaluza-Klein fashion. However, when \t\ falls below 9 3 , the Einstein-dilaton description 
fails and we must employ the fundamental description of strings in order to obtain regular 
behavior at t — 0. 

Let us consider, then, the relevant string action. As we explained in Section IV, for 
winding M2-branes a good gauge choice is A = 1, A 1 = 0. The corresponding gauge-fixed 
worldsheet action was given in JT2J): for p = 2 this reads 

S gf = J drda l -^ v (x»x v - /^tW')) . (30) 

where the fields x^ = (t, x) depend on a and r. This action describes a two dimensional 
field theory with a quartic interaction. 

We confine ourselves to some preliminary remarks about the perturbative behavior of 
()30|) before proceeding to non-perturbative calculations in the next section. In units where 
h is unity, the action must be dimensionless. From the <i-dimensional point of view, the 
coordinates x M have dimensions of inverse mass and /x has dimensions of mass cubed, so a/r 
must have dimensions of mass squared. However, the dimensional analysis relevant to the 
quantization of ()30|) considered as a two dimensional field theory is quite different: the fields 
x^ are dimensionless. From the quartic term we see that fi2@o is a dimensionless coupling, 
suggesting that perturbation theory in /j, 2 6o should be renormalizable. The string tension at 
a time t\ is fi\ = /!2#o|^i|; the usual a' expansion is then an expansion in the Regge slope 
parameter a' = l/(27r/xi), i.e. in negative powers of /i2#o- Conversely, the perturbation 
expansion we are discussing is an expansion in inverse powers of a'. 

These considerations point the way to the resolution of an apparent conflict between 
two facts. Winding M2-brane evolution is as we have seen smooth through t = 0. We 
expect M2-branes to be described by strings near t = 0. The low energy approximation to 
string theory is Einstein-dilaton gravity. Yet, as noted above, in Einstein-dilaton gravity, 
generic perturbations diverge logarithmically with \t\ as t tends to zero. The resolution of 
the paradox is that general relativity is the leading term in an expansion in a' for the string 
theory. As we have shown, however, the good behavior of string theory is only apparent in 
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the opposite expansion, in inverse powers of a'. 

In order to evolve the incoming state defined at large \t\ where general relativity is a 
good description, through small \t\ where string theory is still valid but general relativity 
fails, we must match the standard a' expansion onto the new I /a' expansion we have been 
discussing. We defer a detailed discussion of this fascinating issue to a future publication. 

The action ()30|) describes a string with a time-dependent tension, which goes to zero at t = 
0. There is an extensive literature on the zero-tension limit of string theory (see for example 
(23!), in Minkowski space-time. At zero tension the action is (|30p but without the second 
term. This is the action for an infinite number of massless particles with no interactions. 
Quantum mechanically, there is no central charge and no critical dimension. However, as 
the tension is introduced, the usual central charge and critical dimension appear j^J. 

IX. WORLDSHEET INSTANTON CALCULATION OF LOOP PRODUCTION 

The string theory we are discussing, with action (|HU|) . is nonlinear and therefore difficult to 
solve. We can still make substantial progress on questions of physical interest by employing 
nonperturbative instanton methods. One of the most interesting questions is whether one 
can calculate the quantum production of M2-branes as the universe passes through the big 
crunch/big bang transition. As we shall now show, this is indeed possible through Euclidean 
instanton techniques. 

First, let us reproduce the result obtained in our toy model of winding string production, 
equation (j!8|) . Equation (j!6|) may be re-interpreted as a time-independent Schrodinger 
equation, with t being the coordinate, describing an over-the-barrier wave in a upside- 
down harmonic potential. The Bogoliubov coefficient is then just the ratio of the reflection 
coefficient R to the transmission coefficient T. For large momentum \p\, R is exponentially 
small and T is close to unity since the WKB approximation holds. To compute R we employ 
the following approach, described in the book by Heading j^J. (For related approximation 



schemes, applied to string production, see Refs. |28l l29j.) 



The method is to analytically continue the WKB approximate solutions in the complex 



t-plane. Defining the WKB frequency, u = yP 2 + (/-t#o£) 2 > one observes there are zeros at 
t = ±ip/(fi8o), where the WKB approximation must fail and where the WKB approximate 
solutions possess branch cuts. Heading shows that one can compute the reflection coefficient 
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FIG. 1: The contour for computation of Bogoliubov coefficients in string/membrane production. 

by running the branch cut from t = —ip/(^9o) up the imaginary t axis and out along the 
negative real t axis. Then if one continues an incoming WKB solution, defined below the 
cut, in from t = — oo, below the branch point at t = —ip/([i8o) and back out towards 
t = — oo just above the cut, it becomes the outgoing reflected wave. In the leading WKB 
approximation, the wave is given by w~^e~ l J wdt . Continuing this expression along the 
stated contour, shown in Figure 1, the exponent acquires a real contribution on the parts 
close to the imaginary axis, where t = —it. Hence, the magnitude of the reflection amplitude 
R is given in the first approximation by e^ 2 ^o wdT ^ where the factor of two arises from the 



two contributions on either side of the axis. With w = \jp 2 — (fiO Q r) 2 , one easily sees that 
the exponent agrees precisely with that in (jTHjl . 

A more direct method for getting the exponent in the Bogoliubov coefficient is to start 
not from the field equation (|16|1 but from the original action for the fundamental string, 
f!13|) . We look for an imaginary-time solution (i.e., an instanton) corresponding to the 
WKB continuation described in the previous section. The gauge-fixed particle Hamiltonian 
was given in ((H)) , as H = \ {p 2 — + ((i9 t) 2 ). The corresponding gauge-fixed action is, in 
first order form, 



where as usual dots over a variable denote r derivatives. Before continuing to imaginary 




(31) 



time, it is important to realize that the spatial momentum p = x is conserved (by translation 
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invariance). Hence, all states, and in particular the asymptotic states we want, are labeled 
by p. We are are interested in the transition amplitude for fixed initial and final p, not x, 
and we must use the appropriate action which is not ()31|). but rather 

S = J dr{-p Q i - 'p ■ x - H), (32) 

related by an integration by parts. The p ■ x term contributes only a phase in the Euclidean 
path integral (because p remains real) and the x integration produces a delta function 
for overall momentum conservation. Notice that if we instead had used the naive action 
/ dr^x 2 , we would have obtained a p 2 term in the Euclidean action of the opposite sign. 
Similar considerations have been noted elsewhere j^]. 

Now we continue the action (J3*2|) to imaginary time, setting t = —it E and r = — ir E . 
Eliminating p , the Euclidean action S E = —iS) is found to be 

S E = J dr E Q (i% - {fi9 t E ) 2 + p 2 ) + i'p- xj . (33) 

where dots now denote derivatives with respect to r E . The amplitude we want involves t E 
running from to \p \ and back again: ()33|) is just the action for a simple harmonic oscillator 
and the required instanton is t E = p cos(fi9 T E ) , — n/2 < fi9 r E < n/2. The corresponding 
Euclidean action is 

giving precisely the exponent in pSjl . 

Now we wish to apply this method to calculating the production of winding membrane 
states, described by the action ([3*0J). As in the case of analogous calculations of vacuum 
bubble nucleation within field theory it is plausible that objects with the greatest 
symmetry are produced since non-symmetrical deformations will generally yield a larger 
Euclidean action. Therefore one might guess that the dominant production mechanism is 
the production of circular loops. Let us start by considering this case. The constraint 
n^x 11 ' = implies that the plane of such loops must be perpendicular to their center of 
mass momentum p^. As in the particle production process previously considered, loops 
must be produced in pairs carrying equal and opposite momentum. The Hamiltonian for 
such circular loops is straightforwardly found to be H = | (p 2 + p 2 R — pi + (2iiRfi9ot) 2 ), 
Following the same steps that led to (jSSJ)) we infer that the appropriate Euclidean action is 

S E = J dr E (- (t| + R 2 - (2nn 2 9 Rt E ) 2 + f) + ip ■ . (35) 
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This action describes two degrees of freedom t# and R interacting via a positive potential 
t 2 E R 2 . Up to the trivial symmetries £g — > — tg, R — > —R, there is only one classical solution 
which satisfying the boundary conditions we want, namely starting and ending at = 0, 
and running up to the zero of the WKB frequency function at 2^[iQ^RtE — \p |. This solution 
has Te = R and the Euclidean action is found to be 

my- 
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- / dxVl - x 4 . (36) 

2 JO 



where the last integral is r[|] 2 /(6\/27r), a constant of order unity which we shall denote /. 

The Euclidean action grows like \p\% at large momentum: this means that the total 
production of loops is finite. Neglecting a possible numerical pre-factor in the Bogoliubov 
coefficient, we can estimate the number density of loops produced per unit volume, 

where I is given above. 

From the instanton solution, the characteristic size of the loops and the time when they 
are produced are both of the same order, R ~ \t\ ~ (\p |///2#o)^ ~ (^^o) -1 ^ 3 - The effective 

2 

string tension when they are produced is ~ H2@o\t\ ~ (^2^0)^- 

We have restricted attention so far to the production of circular loops. It is also important 
to ask whether long, irregular strings are also copiously produced. Even though such strings 
would be disfavoured energetically, there is an exponentially large density of available states 
which could in principle compensate. An estimate may be made along the lines of Ref . j3| , 
by simply replacing p 2 with p 2 + fiiN, where \i\ is the effective string tension and N is the 
level number of the string excitations. This picture only makes sense for times greater than 
the string time, so we use the tension at the string time, fx± ~ (^do) 2 ^ 3 - The density of 
string states scales as hence one should replace ()37|) with a sum over N: 

3 - 



Ee^/d d -Ve-(^ +W7V)l/M? ). (38) 



N 

The iV"l beats the \fN so the sum is dominated by modest N, indicating that the production 
of long strings is suppressed. According to this result, the universe emerges at the string 
time with of order one string-scale loop per string-scale volume, i.e. at a density comparable 
to but below the Hagedorn density. 
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Another key question is whether gravitational back-reaction effects are likely to be sig- 
nificant at the transition. As the universe fills with string loops, what is their effect on the 
background geometry? We estimate this as follows. Consider a string loop of radius R in M 
theory frame. Its mass M is 2ttR times the effective string tension H2L, where L is the size 
of the extra dimension. The effective Einstein-frame gravitational coupling (the inverse of 
the coefficient of R/2 in the Lagrangian density) is given by k 2 , = K 2 d+1 / L. The gravitational 
potential produced by such a loop in d spacetime dimensions isj^j 

* = - <w^hs» m 

where A D is the area of the unit -D-sphere, A D = 2nl^±±/r ((D + l)/2). Specializing to the 
case of interest, namely 2-branes in eleven- dimensional M theory, the tension /z 2 is related 
to the eleven dimensional gravitational coupling by a quantization condition relating to the 
four- form flux, reading 

£ = 2n 2 /(n K 2 n ) (40) 

with n an integer. Equations (|39|) and (|40|) then imply that the typical gravitational potential 
around a string loop is 

1 fl^i 

• = -64^—^"'-^" (41) 

up to numerical factors. 

We conclude that the gravitational potential on the scale of the loops is of order 9q and 
therefore is consistently small for small collision rapidity. Since the mean separation of 
the loops when they are produced is of order their size R, this potential $ is the typical 
gravitational potential throughout space. Multiplying the tt component of the background 
metric (J2J) by 1 + 2$ and redefining t, we conclude that the outgoing metric has an expansion 
rapidity of order ~ 8q(1 + C6q) with C a constant of order unity. We conclude that for small 
#0 the gravitational back-reaction due to string loop production is small. Note that loop 
production is a quantum mechanical effect taking place smoothly over a time scale of order 
(//2#o)~ 3 - Therefore if the rapidity of the outgoing branes alters as we have estimated, it 
happens smoothly and not like the jump in #0 discussed at the end of section VI. Therefore 
the picture of loop production is consistent with the comments made there. 



26 



X. STRINGS ON LORENTZIAN ORBIFOLDS ARE NOT REGULAR 



We have shown that strings constructed as winding M2-branes on M.q x R 9 are analytic 
in the neighborhood of t = 0. This was the setup originally envisaged in the ekpyrotic model, 
where collapse of the M theory dimension was considered. Subsequently, a number of authors 
investigated the simpler case of string theory on A4c x R 8 , considered as a Lorentzian orbifold 
solution of ten dimensional string theory. This is a simpler, but different setting, hence we 
expressed misgivings Q,ll8| about drawing conclusions from these reduced models. We shall 
now explain why the behavior in the Lorentzian orbifold models is significantly worse than 
in M theory and, hence, why no negative conclusion should be drawn on the basis of the 
failed perturbative calculations. 

Consider string theory on the background (J2J). Let us choose the gauge x° = t = r, and 
g^x^x'^ = 0. In this gauge we express the string solution as 9(t, a) and x(t, a). If a classical 
solution in this gauge possesses a singularity at t = in the complex t-plane, for all a, then 
there can be no choice of worldsheet coordinates r and o which can render the solution 
analytic in the neighborhood of t = 0. For if such a choice existed, one could re-express r in 
terms of t and, hence, x(t, a) would be analytic. We shall show that generically, for strings 
on Lorentzian orbifolds, the solutions possess logarithmic singularities i.e., branch points, 
rendering them ambiguous as one circumvents the singularity in the complex t-plane. 

The demonstration is straightforward, and our argument is similar to that in earlier 
papersQ]. We are only interested in the classical equations of motion and we may compute 
the relevant Hamiltonian from the Nambu action, 



H = J doH = J do 



n 2 7T 2 



\ — 2 +- + ^{{t6'Y + x'i). (42) 



The Hamiltonian equations allow generic solutions in which 7Tq(o) tends to any function of a 
as t tends to zero. From its definition, the Hamiltonian density TC then diverges as t . The 
equation of motion for 9 is 9 = ng/(t 2 T-C), implying that 9 — > ±t _1 independent of a. This 
implies a leading term 9 ~ ±logt, independent of a. Recalling that M.c may be rewritten 
in flat coordinates by setting T = t cosh (9 and Y = tsinh^, one readily understands this 
behavior. A geodesic in the (T, Y) coordinates is just Y = VT, with V a constant. At small 
t this requires e e or e~° to diverge as which is just the result we found. 

We conclude that generic solutions to the Hamiltonian equations possess branch points at 
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t = meaning that the solutions to the classical equations are ambiguous as one continues 
around t = in the complex t-plane. This is a much worse situation that encountered for 
winding branes in M theory. 

The second problem occurs when we quantize and construct the associated field theory. 
Then the bad behavior at t = corresponds to a diverging energy density which renders 
perturbation theory invalid. Since the previous problem can be seen even in pointlike states, 
let us focus attention on those. As we discuss in detail in Appendix 4, one needs to use 
the metric on the space of coordinates in order to construct the quantum Hamiltonian. In 
this case, the metric on the space of coordinates is the background metric, (j2J. The field 
equation for point particles is then given from (J66j) : Fourier transforming with respect to x 
and 9, it reads 

4> + U=-P<f>-%, (43) 

n 

where kg is quantized in the usual way. This is the equation studied in earlier work[3| on 
quantum field theory on M.q x -R^™ 1 - The generic solutions of ()43|) behave for small t as 
log t for kg = 0, or t ike for kg ^ 0. In both cases, the kinetic energy density <fi 2 diverges as 
t~ 2 . Similar behavior is found for linearized vector and tensor fields on M. c x R d ~ l . These 
divergences lead to the breakdown of perturbation theory in classical perturbative gravity, 
an effect which is plausibly the root cause of the bad behavior of the associated string theory 



scattering amplitudes 



_Dlv 



As we have stressed, in the sector of M theory considered as a 
theory of membranes, describing perturbative gravity, this effect does not occur. The field 
equations are regular in the neighbourhood of t = and there is no associated divergence 
in the energy density. 

With hindsight, one can now see that directly constructing string theory on Lorentzian 
orbifolds sheds little light on the M theory case of interest in the ekpyrotic and cyclic models. 
Whilst the orbifolding construction provides a global map between incoming and outgoing 
free fields, it does not avoid the blueshifting effect which such fields generically suffer as they 
approach t — 0, which seems to lead to singular behavior in the interactions (although Ref. 
14j argues that a resummation may cure this problem). 

An alternate approach involving analytic continuation around t — has been simul- 
taneously developed, but so far only implemented successfully in linearized cosmological 
perturbation theory P. 11^. This method may in fact turn out to work even in the M theory 
context. The point is that by circumnavigating t = in the complex t-plane, maintaining 
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a sufficient distance from the singularity, one may still retain the validity of linear pertur- 
bation theory and the use of the Einstein equations all the way along the complex time 
contour. The principles behind this would be similar to those familiar in the context of 
WKB matching via analytic continuation. 

In any case, the main point we wish to make is that now that we have what seems like a 
consistent microscopic theory for perturbative gravity, valid all the way through t = 0, we 
have a reliable foundation for such investigations. 

XI. CONCLUSIONS AND COMMENTS 

We have herein proposed an M theoretic model for the passage through a cosmological 
singularity in terms of a collision of orbifold planes in a compactified Milne Aic x R 9 
background. The model begins with two empty, flat, parallel branes a string length apart 
approaching one another at constant rapidity 8q. With this initial condition, we have argued 
that the excitations naturally bifurcate into light winding M2-brane modes, and a set of 
massive modes including the Kaluza-Klein massive states. It is plausible that the massive 
modes decouple as their mass diverges. The light modes incorporate perturbative gravity 
and, hence, describe the space-time throughout the transition. Our finding that they are 
produced with finite density following dynamical equations that propagate smoothly through 
the transition supports the idea that this M-theory picture is well-behaved and predictive. 
Our model also suggests a string theory explanation of what goes wrong with Einstein 
gravity near the singularity: Einstein gravity is the leading approximation in an expansion 
in a', but the winding mode picture is a perturbative expansion in I /a'. 

Our considerations have been almost entirely classical or semi-classical, although we 
believe the canonical approach we have adopted is a good starting point for a full quantum 
theory. Much remains to be done to fully establish the consistency of the picture we are 
advocating. In particular, we need to understand the sigma model in 1)46)1 . and make sure 
that it is consistent quantum mechanically in the critical dimension. Second, we need to 
learn how to match the standard a' expansion to the I /a' expansion which we have argued 
should be smooth around t = 0. Third, we need to incorporate string interactions. Although 
the vanishing of the string theory coupling around t = suggests that scattering plays a 
minor role, the task of fully constructing string perturbation theory in this background 
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remains. 

Assuming these nontrivial tasks can be completed, can we say something about the 
significance of this example? The most obvious application is to the cyclic and ekpyrotic 
universe models which motivated these investigations, and which produce precisely the initial 
conditions required. As has been argued elsewhere Q, Q], the vicinity of the collision is 
well-modeled by compactified Milne times flat space. Within this model, the calculations 
reported here yield estimates of the reheat temperature immediately after the brane collision 
i.e. at the beginning of the hot big bang. Furthermore, if our arguments of Section VIII are 
correct, they should in principle provide complete matching rules for evolving cosmo logical 
perturbations through singularities of the type occurring in cyclic/ekpyrotic models. As we 
discussed briefly in Section X, it is plausible that this matching rule will recover the results 
obtained earlier for long-wavelength modes by the analytic continuation method 
although that remains to be demonstrated in detail. 

Many other questions are raised by our work. Can we extend the treatment to other 
time-dependent singularities? The background considered here, M. c x R 9 is certainly very 
special being locally flat. Although the string frame metric is singular, it is conformally flat. 
One would like to study more generic string backgrounds corresponding to black holes, or 
Kasner /mixmaster spacetimes in general relativity. As we have argued, there is no reason to 
take the latter solutions seriously within M theory since they are only solutions of the low 
energy effective theory, which fails in the relevant regime. Nevertheless, they presumably 
have counterparts in M theory, and it remains a challenge to find them. We would like 
to believe that by constructing one consistent example, namely M theory on Aic x R 9 we 
would be opening the door to an attack on the generic case. 

Such a program is admittedly ambitious. However, should it succeed, we believe it would 
(and should) completely change our view of cosmology. If the best theories of gravity allow 
for a smooth passage through time-dependent singularities, this must profoundly alter our 
interpretation of the big bang, and of the major conceptual problems of the standard hot 
big bang cosmology. 
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Appendix 1: Canonical Treatment of p-Branes in Curved Space 

In this Appendix we review the canonical treatment of p-Brane dynamics in curved space. 
A similar approach is taken in the recent work of Capovilla et al. j^j. Earlier treatments 



include Refs. 
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We start from the action in Polyakov form. No square roots appear and the ensuing gen- 
eral Hamiltonian involves only polynomial interactions. The action for a p-brane embedded 
in a background space-time with coordinates x^ and metric g^ v is: 

S = ~p p J d p+1 a^ (l^d^d^g^ - (p - 1)) , (44) 

where p, v is a mass per unit p- volume. The p-brane worldvolume has coordinates a a , a = 
0,1, ... ,p in which the metric is 7 aj g. 

One way to proceed is to vary (j44j) with respect to 7^, hence, obtaining the constraint 
expressing the worldvolume metric j a p in terms of the the induced metric d a x^dpx v g^ u . 
Substituting back into (jH| one obtains the Nambu action for the embedding coordinates 
x^{a) i.e., —p p times the induced p-brane world volume. 

One can do better, however, by not eliminating the worldvolume metric so soon. Instead it 
is better to retain the as independent variables and derive the corresponding Hamiltonian 
and constraints corresponding to evolution in worldvolume time r. It is convenient to express 
7 a/ 3 in the form frequently used in canonical general relativity: the worldsheet line element 
is written 

"j a pda a da p = (-a 2 + p k p k )dr 2 + l^rdo 1 + T^-doW, (45) 

where a 1 , i = 1, . . . ,p are the spatial worldvolume coordinates, (3 k is the shift vector and a is 
the lapse function. The good property of this representation is that the metric determinant 
simplifies: 7 = — a 27 f. 

In the following discussion we shall for the most part assume p > 1 and then comment 
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on the amendments needed for p = 0,1. Using (|4*5j) the action ()44|) becomes 

If i 

S = — J drcPaarfz 



- 2 x^g, u + 2^x^ + (f j - ^Xifig,., - (p - 1) 



(46) 



The canonical momenta conjugate to x M are found to be 



i 

7T,i = H P —{x v ~ P l x v ,d9mv (47) 

No time derivatives of a, (5 l and 7^ appear in the Lagrangian, hence, the corresponding 
conjugate momenta vanish. In Dirac's language Q], these are the primary constraints. 

7l a PS 0, 7Tj Rj 0, 7Ty ~ 0. (48) 

Since Poisson brackets between momenta vanish, these constraints are first class. The total 
Hamiltonian H then consists of the usual expression H = J d p a [tiax a — L] , 



H = d p a 



a xTT^cT + /fc#r„ + (r j <iXj9^ - (P - 1)) 



(49) 



2^72 " " 2 

plus an arbitrary linear combination of the primary first class constraints 

Additional secondary constraints are obtained from the Hamiltonian equations for 7t a , -k, 
and TTij. Insisting these vanish as they must for consistency with (}4*%|) . one finds 



C = + //Jy (fM 2 ^ - (p - 1)) « 0, 

Cj = x^ir^ ~ 0, 

Cij = 7ij - x ^i x j9^ ~ 0. (50) 

Following Dirac's procedure, we can now try to eliminate variables using second class con- 
straints (constraints whose Poisson brackets with the other constraints does not vanish). 
In particular, it makes sense to eliminate 7^ since the corresponding momenta vanish 
weakly. It is easy to see that CV,- are second class since their Poisson brackets with 717,- are 
nonzero. Hence, we eliminate 7^ using the constraint, obtaining 

C ee ^ v tT + /ijDet(x%.^) ~ 0, Cj ee xfa ^ 0, (51) 

as our new constraints, on the remaining variables x M and 717. The matrix x^x^g^ is the 
induced spatial metric on the brane. When written this way, the C and Ci constraints 
have zero Poisson brackets with the remaining primary constraints 7i a and 77. A lengthy 
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but straightforward calculation establishes that all Poisson brackets between C and Cj are 
weakly vanishing (see Appendix 2) and, hence, that we have a complete set of first class 
constraints consisting of n ai m, C{ and C. 

The canonical Hamiltonian (J49|) is now seen to be a linear combination of C and Cj, with 
coefficients depending on a and f3 l respectively. The general Hamiltonian consists of a sum of 
this term plus an arbitrary linear combination of the first class constraints, / da(v a 7r a + v l 7Ti) 
where v a and V{ are arbitrary functions of the worldvolume coordinates. From Hamilton's 
equations, one infers that a = v a , and f3 l = v l . Therefore, a and (3 are completely arbitrary 
functions of time. As Dirac emphasizes, one can then forget about a, 7r a , f3 t and 7Tj and just 
write the total Hamiltonian for the surviving coordinates as x M and 7r M as 



i.e. a linear combination of the constraints (|51j) with arbitrary coefficients A and A 1 . Dif- 
ferent choices of A and A 1 then correspond to different choices of worldvolume coordinates. 

We count the surviving physical degrees of freedom as follows. We start with the 2(d + 1) 
coordinates and momenta tt 11 , each functions of the p + 1 worldvolume coordinates. Then 
we impose the p+1 constraints C = Cj = 0. Finally, in order to specify time evolution we 
must pick p+1 arbitrary functions A and A{. The remaining physical degrees of freedom 
are 2(d + 1) — 2(p + 1) = 2{d — p) in number, the right number of transverse coordinates 
and momenta for the p-brane. 

We have tacitly assumed p > 1 in the above analysis. The following minor amendments 
are needed for p = and 1. For p = 0, the worldvolume metric involves a only and one 
can ignore anything with an i index except the determinant, which is replaced by unity. In 
particular there is no integration over a in the Hamiltonian, and the canonical momentum 
density 7r M is replaced by the momentum p^. The only constraint is 



which is just the usual mass shell condition. The general Hamiltonian consists of an arbitrary 
function of r times C. 

For p — 1 the action is invariant under conformal transformations of the worldvolume 
metric, and hence only two independent combinations of the three worldvolume metric 
variables appear in the decomposition ([45)1 . The corresponding two momenta vanish and 




(52) 



C= Ptl p u g^ + p, 2 ^0, 



(53) 
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these are the primary constraints. Through Hamilton's equations one finds the following 
secondary constraints: 



C = w v gT + filxVg^ « 0, 

d EE TT^' W 0, (54) 

where primes denote derivatives with respect to a 1 ee a. The general Hamiltonian again 
takes the form (|52|) . with p = 1. 

Appendix 2: Poisson bracket algebra of the constraints 



In the canonical theory , one considers arbitrary functions of the canonical variables, 
and of the time r. In our case, the canonical variables are fields x^{a) and ^(u) depending 
upon a, which is regarded as a continuous index labeling an infinite number of canonical 
variables. In particular, the constraints in (j51|) are infinite in number. In this Appendix we 
shall show that the Poisson bracket algebra of the constraints closes, and, hence, that, in 
Dirac's terminology, they are first class. 

The Poisson bracket between any two quantities M and N, which may be arbitrary 



functions of the canonical variables (local or nonlocal in a) and of the time r, is given by 

f m f 9M ON ON dM \ 
{M,H) S **{ - , (55) 



where (dx^(a')/ dx u (a)) = (d7r u (a')/d7r tl (a)) = 5%5 p (a — &'), with other partial derivatives 
being zero. 

One way to calculate the Poisson brackets between a set of constraints C and Cj, is to 
start from a putative Hamiltonian 

H = J dPai^C + A'a), (56) 

where A and A % are arbitrary functions of a, and then compute Hamilton's equations for the 
r derivatives of x M and tt^. We then use these to determine the corresponding r derivatives 
of C and Cj. Setting these equal to {C, H} and {Cj, H} with H given by (fHHjl. we are able 
to infer the Poisson brackets between the constraints. For the Hamiltonian ()56|) with C and 
Ci given in (|51|) . Hamilton's equations read 

x^ = Ag" u 7r u + 
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(57) 



where dots denote r derivatives and 7^ = x^x^g^ is the induced spatial metric on the 
brane and 7 its determinant. We have made use of the formula dry = J^dry^. 

Using (|H7j) , it is a matter of straightforward algebra to compute C and C; and hence infer 



all of the Poisson brackets. We find 



{C(a),CM')} 



{C{a),C{a')} 




{Ci(a),Cj(a')} 



(58) 



The right hand side consists of linear combinations of the constraints and, hence, it vanishes 
weakly. We conclude that the constraint algebra closes and, hence that the constraints are 
first class. Notice that the case of strings, p — 1, is specially simple since 7 7 11 = 1 and the 
Poisson bracket algebra is linear, with field-independent structure constants. 

The calculation also provides a consistency check on our Hamiltonian ©, which is pre- 
cisely of the form (|56p. since it implies the constraints are preserved under time evolution 
in r. 

Appendix 3: Equivalence of gauge-fixed Hamiltonian and Lagrangian equations 

In this Appendix we establish that the Lagrangian equations following from the gauge- 
fixed action ()30)1 for winding M2-branes are equivalent to the Lagrangian equations for 
a string in the time-dependent background g^ u = \t\r]^ u , in a certain string worldsheet 
coordinate system. This is in accord with our general arguments. 

The equations of motion following from the gauge-fixed action (j3*Uj) are: 



i = tx' 2 + tt' 2 + t"t 2 . 



(59) 



and the constraints take the form 





(60) 
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We want to compare these equations with the Lagrangian equations of motion following 
from the Polyakov action ©, with p — 1. These are 

+{-l) l2 Kx{l TT drX v d T x x + r a d a x v d a x x ) = 0, (61) 

where T^ x is the Christoffel symbol for the background metric. 

We also have the constraints that the worldsheet metric 7^ is conformal to the induced 
metric on the string. We have the freedom to choose worldsheet coordinates on the string, 
but since the equations are conformally invariant, only the conformal class matters. The 
choice 7 Q/ g = f2 2 diag(— t 2 , 1) is found to yield the two constraints (Jot^ . 

For our background, g^ u = [tlr]^, we have nonzero Christoffel symbols r{] = l/(2t), 
r* Q = T l j = 5ij/(2t), r°. = 5ij/(2t), where i runs over the background spatial indices 1 to 
d — 1. The string equations of motion (J61)) are then found to be equivalent to (J59|) . for all 
nonzero t. 

From the string point of view, this choice of gauge would seem arbitrary, and indeed it 
would appear to be degenerate at t — 0. Yet, as we have seen, this gauge choice is just 
A = 1 and A 1 = 0, which is entirely natural from the canonical point of view. It has the 
desirable property that the equations of motion and the constraints are regular at t — 0, and 
from the general properties of the canonical formalism we are guaranteed the existence of 
an infinite class of coordinate systems, related by nonsingular coordinate transformations, 
in which the equations of motion will remain regular. 

Appendix 4: Ordering ambiguities and their resolution for relativistic particles 

In the main text we have discussed the canonical Hamiltonian treatment of relativistic 
particles and p-branes. When one comes to quantize these theories in a general back- 
ground, certain ordering ambiguities appear which must be resolved. Here we provide a 
brief overview, following the more comprehensive discussion in Ref. 

The field equation for a relativistic particle is simply the expression of the quantum 
Hamiltonian constraint H = 0, in a coordinate space representation. The first task is to 
determine the representation of the momentum operator in this representation, and then 
that of the Hamiltonian operator H. As we shall now discuss, this requires knowledge of 
the metric on the space of coordinates. We shall only deal with the point particle case. 
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The classical Hamiltonian constraint for a massive particle in a background metric g^ u 
reads 

g^Ptfflv + m<2 ~ 0. (62) 

First we attempt to determine the coordinate space representation of p^, consistent with the 
quantum bracket: 

\x^p v \=iM». (63) 

One choice is p^ = —ifid^ but this is not unique: the representation p^ = —ih(d^ + / M ), with 
fn any function of the coordinates x M and r, is equally good as far as (jSHj) is concerned. 

We now show how / M may be determined from the additional requirement that p^ be 
hermitian i.e., that the momentum be real. In the coordinate space representation, this 
requirement reads 

(X\PM = I d d x{-g{x))^ ( X > M 0) = (<j>\p„\xY = J d d x(-g(x))i (<f>* PflX Y (64) 

where the integration runs over the space of coordinates and g^ v is the metric on that space. 
It is straightforward to check that the naive operator —ihd^ is in fact not hermitian for 
general g^, but that 

?„ = -a (a„ + \ (a>(-s))) = -in(- g )-ia,(- g )i (65) 

is. This discussion uni q oely determines the re* part of /„: - imaginary part may be 
absorbed in an unobservable phase of the wavefunction[15j. 

Similarly, when we consider the Hamiltonian constraint (|62|). the questions arise of where 
to place the g^ u relative to the p^s, and whether to include any factors of the metric deter- 
minant g. The resolution is familiar: if we write the Hamiltonian in covariant derivatives on 
the space of coordinates, it will be hermitian since we can integrate by parts ignoring the 
\J—g factor in the measure. This suggests setting the first term in (}6*2*|) equal to the scalar 
Laplacian: 

sTPvPu - -n\-g)^d,{-g)^g^d, = (-gylp^-g) 1 * g^ Pv (- g yl . (66) 

It is straightforward to check that this is the only choice of ordering and powers of (—g) 
which is hermitian and has the correct classical limit. Nevertheless, this ordering is not 
immediately apparent! More generally, one can also include terms involving commutators of 
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p M which are zero in the classical limit, but which produce the Ricci scalar R in the quantum 
HamiltonianQj. In the space-time we consider R is zero. Hence, such terms do not arise. 
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